Bekki-Nozaki Amplitude Holes in Hydrothermal Nonlinear Waves 
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We present and analyze experimental results on the dynamics of hydrothermal waves occuring in 
a laterally-heated fluid layer. We argue that the large-scale modulations of the waves are governed 
by a one-dimensional complex Ginzburg-Landau equation (CGLE). We determine quantitatively 
all the coefficients of this amplitude equation using the localized amplitude holes observed in the 
experiment, which we show to be well described as Bekki-Nozaki hole solutions of the CGLE. 
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The status and nature of the so-called amplitude equa- 
tions which can be derived in the vicinity of symmetry- 
breaking instabilities is now well-established They 
are "universal" in so far as they essentially depend on 
the symmetries of the physical system and of its bifur- 
cated solutions, but also because they often remain valid, 
at least at a qualitative level, even far away from the 
instability threshold However, determining accu- 

rately the coefficients of the underlying relevant ampli- 
tude equation from experimental data remains a difficult 
task, especially in these far-from-thrcshold regimes. 

The complex Ginzburg-Landau equation (CGLE), 
which describes the large-scale modulations of the bi- 
furcated solutions near oscillatory instabilities, is per- 
haps the most-studied amplitude equation |^ . This priv- 
iledged situation is due to both its relevance to many ex- 
perimental situations and to the variety of its dynamical 
behavior, in particular its spatiotemporal chaos regimes. 
One of the landmarks of the CGLE is that it possesses 
localized "defect" solutions. Even in one space dimen- 
sion, where no topological constraint exists, numerical 
simulations of the CGLE and analytical work 
have revealed the existence and importance of various 
amplitude hole solutions, which can often be seen as the 
"building blocks" of the complex spatiotemporal dynam- 
ics observed. In particular, the one-parameter family of 
traveling hole solutions discovered by Bekki and Nozaki 
1^ has been shown to play an important dynamical role 
in a large portion of parameter space including in regions 
where they are linearly unstable j^. Similar objects have 
been identified in various experimental contexts of a pri- 
ori relevance, e.g. Rayleigh-Benard convection and cou- 
pled wakes However, to our knowledge, there is still 
no case where a direct comparison with known solutions 
of the CGLE could be achieved. 

In this Letter, we present a quantitative comparison 
of localized amplitude holes observed in an experiment 
with hole solutions of the CGLE, the relevant amplitude 
equation. We use the observed holes to fully determine 
the coefficients of the underlying CGLE. This provides 
clear-cut evidence of Bekki-Nozaki holes in an experi- 
mental context. Our system is a long, straight, and nar- 
row convection cell in which a thin fluid layer with a free 
surface is subjected to a horizontal temperature gradi- 
ent. Hydrothermal nonlinear waves appear via a direct 



Hopf bifurcation, indicating the relevance of the CGLE. 
The spatiotemporal dynamics of the waves exhibits local- 
ized amplitude holes. The basic scales of the equivalent 
CGLE are determined using the regular part of the wave 
trains. Data collected in the vicinity of amplitude holes 
show that they have the structure of Bekki-Nozaki so- 
lutions. This also provides estimates of the remaining 
coefficients of the CGLE, an approach which, we argue, 
could be efficient in other experimental contexts. Finally, 
the overall consistency of our results is checked. 



T+ 


Air 

> 


T_ 


Fluid 



-i 











T"! , ■ ■ 1 ■ 




source ^ 



FIG. 1. Experimental cell and basic wave pattern. Above: 
schematic side view. Below: top view and instantaneous shad- 
owgraphic image of the wave pattern. 

The experimental setup is schematically described in 
Fig. 0. A layer of fluid (silicon oil of viscosity v = 0.65cSt 
and Prandtl number P = 10) of height h is confined be- 
tween two copper blocks maintained at fixed tempera- 
tures Tjf- and r_ by thermostated water circulation, and 
a bottom glass plate. This forms a straight, narrow chan- 
nel of length Lj. = 25cm and width Ly = 2cm. As soon 
as the temperature difference AT = T-|_ — T_ is not zero, 
a basic flow sets in. It consists of a surface flow towards 
the cold side with a bottom recirculation. Increasing AT, 
the basic flow becomes unstable to traveling hydrother- 
mal waves |l^] via a supercritical Hopf bifurcation . 
We observe these waves by low-contrast shadowgraphy, 
which captures the vertical average of the temperature 
gradient variations (surface waves exist, but their effect is 
negligible). In this geometry, the waves propagate away 
from a "source region" located arbitrarily on the cold 
wall, and the end boundaries at a; = 0, Lx act as sinks 
with no apparent reflection (Fig. For h = 1.2mm, cor- 
responding to the experiment reported below, the source 
region emits curved waves which become planar further 
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away (Fig. bottom) and propagate along the x-axis. 
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FIG. 2. Spatiotemporal evolution of the wave pattern at 
h = 1.2mm and AT = 5. IK (yielding a Marangoni number 
Ma ~ 950). Only a central portion of length 18.5cm is shown 
during 100s. (a): Original data, (b): Evolution of the mod- 
ulus \A\ (black: \A\ = 0; white: \A\ = 1) and (c): phase 
gradient k (dark: fc < 0; bright: k > 0). 

Figure ^ja presents a typical spatiotemporal evolution 
as obtained from the acquisition, with a fixed-gain cam- 
era, of a single 512-pixel line (of negligible width) along 
the X-axis in the center of the cell. Here, the source 
appears as a rather ill-defined, erratic object. (Closer 
to the Hopf bifurcation, steady, regular, evolution is ob- 
served.) Fourier analysis of diagrams such as Fig. ||a 
reveals that on each side of the source only waves prop- 
agating away from the source are present and that they 
are approximately monochromatic (the second harmonic 
is two orders of magnitude smaller). More precisely, re- 
stricting ourselves to one side of the source (say x > 90 
on Fig. 1^), we can write the recorded physical variable: 

V{x,t) — A{x,t) eyip[i{kox — uJot)] + c.c. (1) 

where fcg is the dominant wavelength, loq the basic fre- 
quency, and A a one-dimensional complex field describing 
the (large-scale) modulations of this wave. (On the other 
side of the source, one has to change the sign of fcp.) 

Using complex demodulation techniques, A can be ex- 
tracted from the experimental data. Figs. ||b,c show 



the spacetime evolution of \A\ and k = g ± fco where 
q — dj:aig{A). In these pictures, the localized deforma- 
tions of the waves visible in Fig. ^ clearly appear as 
propagating amplitude holes across which the phase gra- 
dient varies rapidly. At some space-time points, \A\ even 
vanishes and the phase gradient diverges: a space-time 
dislocation occurs (Fig. 0a). The amplitude holes can 
be seen as the objects mediating the evolution to wave 
patterns more regular than those emitted by the source. 

Our system clearly calls for a one-dimensional model. 
The waves arise via a supercritical Hopf bifurcation. 
Away from the source, they propagate only in one di- 
rection. All this indicate that the evolution of A could 
be governed by a single CGLE on each side of the source, 
even though the regime studied here takes place at finite 
distance from threshold (for h = 1.2mm, ATc ~ 4.3K, 
and thus the relative distance to threshold is £ = 0.19 for 
AT = 5. IK). We thus suppose that A obeys: 

To{dt + Virdx)A = 

eA + (^{l+ia)d^,A - g{l+iP)\A\^A (2) 

where Vg is the group velocity of the waves, tq and 
are the basic time- and length-scales of the wave modu- 
lations, and (7 is a real number. Below, we estimate, from 
the data of Fig. ^, all the coefficients of Eq. (H) and check 
the overall consistency of our hypothesis. 

The linear part of the variation of the local frequency 
uj with the local wavenumber k yields our estimate of 
the group velocity: Vg = div/dk ~ — 1.16mm/s (Fig. |^). 
This is consistent with the average value of the velocity 
of small perturbations, estimated at —1.15 ± 0.25mm/s, 
to be compared to the phase velocity = 1-^0/^0 — 
— 2.8mm/s. This confirms that the source is indeed a 
source, since perturbations do propagate outward. 

Fig. ^ shows the variation of \A\ with k as deter- 
mined from the portion of Fig. ^ at the left of the source 
(x < 90mm). The maximum amplitude is observed for 
the basic wavenumber: kg ~ — l.llmm^^. Space-time 
points away from the localized amplitude holes corre- 
spond to the large \A\ (say \A\ > 0.5) portion of the 
curve. Locally around these points, the solution of (||) is 
expected to be close to one of the phase- winding solutions 
of wavevector q = k + ko (see e.g. Q): 

A = Aq exp[i(gx - Lv^t)] with A^ ^ {e - ilq^)/g 
and ujq = [el3+{a- l3)ilq^]/To~Vgq (3) 

The linear variation of |Ap with q^ is confirmed in 
Fig. Ic, yielding ^o/Ve 2.53mm and ./e/g ~ 0.00054 
(a.u.). Note that we thus have 3> £.o/V^ ~ 3/fco: the 
cell is effectively "infinite" and the variations of A occur 
on scales significantly larger than the basic length k^^. 

Timescale tq can be estimated from the real part of 
the spatial linear growth rate of waves near the source, 
which is equal to e/(roWg) [Q. From Fig. we find 
To/e ~ 8.5 ± 0.5s, about four times the basic period 
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2tt/u!q ~ 2.03s, confirming that the variations of A are 
slow compared to the basic oscillations. Note that tq is of 
the same order as the viscous diffusion time /v — 2.2s. 
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FIG. 3. Analysis of the data of Fig. ^ situated at the left 
of the source {x < 90mm). (a-c): for each bin around a 
given k value, all the corresponding space-time points were 
determined, and averages calculated on each bin. (a): local 
frequency u; vs local wavenumber k; the slope of a linear fit 
gives the group velocity Vg. (b): local amplitude \A\ vs k; 
solid line: see text, (c): \A\'^ vs (same data as in (b)). (d): 
two instantaneous profiles of |j4|(a;) near the source taken at 
time and 30s; insert: In \A\ vs x and linear fit (thin line). 

At this stage, all the basic scales of Eq. (||) have been 
estimated. To determine the remaining two parameters 
a and /3, global quantities deriving from the "wave part" 
of the data could, in principle, be sufficient. For example. 
Fig. ||a could be used to extract the expected variation of 
LUq with q. But the data is too noisy to yield any meaning- 
ful estimate of a and /3. Moreover, as long as the source 
is not controlled, the "input" waves cannot be varied at 
will to explore the family of solutions (||), contrary to 
other experimental situations We now focus instead 
on the localized amplitude holes already mentioned. 

Many localized, propagating objects connecting two 
phase-winding solutions have been observed numerically 
in the one-dimensional CGLE Analytical methods 
are largely limited, so far, to solutions depending only 
on the reduced variable ^ — x — Vht where Uh is the (con- 
stant) velocity of the object 1^,1). Using this Ansatz, the 
CGLE reduces to a third-order ordinary differential equa- 
tion (ODE) whose fixed points are the phase-winding so- 
lutions (^. Localized objects connecting two such solu- 
tions of wavevector and qn appear, within the Ansatz, 
as homoclinic (^l = ^r) or heteroclinic ((/l qn) orbits. 

The holes observed in Fig. |2| are not stable structures 
connecting two infinite phase winding solutions, but they 



subsist long enough, and can be sufficiently isolated to re- 
veal that their wings are indeed well described by phase 
winding solutions. (As a matter of fact, we repeatedly 
assimilated, above, the large-amplitude regions separat- 
ing the holes to portions of these solutions.) Fig. ^ shows 
an isolated hole extracted from Fig. |[ One can mea- 
sure rather accurately the two wavenumbers q^ — 0.13 
and qb. = —0.32 (in the CGLE frame) connected by the 
central hole, which can thus be tentatively seen as an 
heteroclinic orbit in the ODE Ansatz. 

All sufficiently localized structures on Fig. ^ also con- 
nect two different wavenumbers. This rules out the ho- 
moclinic holes recently studied by van Hecke and 
leave, as possible candidates, the family of hole solutions 
found by Bekki and Nozaki . The explicit form of these 
solutions is too lengthy to be given here (see, e.g., Q). 
They form a one-parameter family (at fixed a, j3) which 
can be parametrized by, e.g., the velocity Uh of the hole. 
They take the shape of an exponentially- localized ampli- 
tude hole with a minimum amplitude | accompanied 
by a rapid phase shift a. 
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FIG. 4. Amplitude hole. Solid lines: experimental data 
extracted from a cut of Fig. |^ at f = 60s. Dashed lines: corre- 
sponding Bekki-Nozaki hole solution of the equivalent CGLE. 
(rescaled variable x' = xv^/^o). 

To compare data such as that of Fig. ^ to these solu- 
tions, we need to determine the values of a and (3 and 
the "optimal" solution of the corresponding family. We 
proceed as follows: we estimate gL, 9r, and |v4|niin from 
the data since we found these were the characteristics of 
the hole for which the most accurate measurement can 
be made. We find, for all values of the (a, (3) plane where 
it exists, the Bekki-Nozaki hole solution with the mea- 
sured value of gL + <1k- We then select the (codimension 
1) subsets of the (a, (3) plane where, moreover, this hole 
solution possesses the measured value of or the esti- 
mated value of |A|i„in (Fig. ||, dashed lines). These two 
lines intersect, yielding the desired values of a and /3. 
Taking into account the error bars on gL, 5r, and |A|min, 
we find a ^ -1.5 ± 0.5 and (3 = -0.4 ± 0.05. By the 
same token, the hole solution is also uniquely determined 
(Fig. H, dashed hues). Its velocity Wh, its width, its phase 
shift a are all found consistent with the data. 

We repeated the procedure for other amplitude holes 
found in Fig. ^ We found almost the same a and P val- 
ues, to the accuracy of the estimates. This strengthens 
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the confidence in the results, since different objects mov- 
ing at different velocities, connecting different wavenum- 
bers all yield the same parameter values. We also per- 
formed a final check, by plotting, for the estimated values 
of a and /?, the variation of |A|inin with the local phase 
gradient at the "bottom" of the hole along the family of 
solutions (Fig. ||b, solid line). The agreement with the 
small- 1^1 values measured on Fig. ^jis very good. This is 
an additional indication that all the low-amplitude points 
are indeed located "inside" Bekki-Nozaki holes. 
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FIG. 5. Parameter plane of tlie CGLE. STl: spatiotempo- 
ral intermittency regions (from [^). B-F: Benjamin- Feir line 
below which all solutions (^) are unstable. Dashed lines: see 
text (grey areas: error bars). Circle: estimated {a, (3) values. 

For the estimated values of a and (3, the CGLE is in the 
parameter region where the phase-winding solutions (|^) 
are linearly stable for \q\ small enough and no sustained 
spatiotemporal disorder exists in one space dimension . 
Moreover, the Bekki-Nozaki amplitude holes are linearly 
unstable This is not in contradiction with the dy- 

namics observed in the experiment: the Bekki-Nozaki 
holes, although unstable, exist, and can constitute im- 
portant building blocks of even chaotic dynamics ||] . The 
waves emitted by the source can be locally attracted to 
this family of unstable fixed points before escaping along 
its unstable manifold (a mechanism also invoked by van 
Hecke in Q ) . The tendency of the waves trains to become 
more regular away from the source (see Fig. ||) is consis- 
tent with disorder being only transient in the CGLE with 
the estimated parameter values. 

In summary, we have presented experimental results on 
the dynamics of the nonlinear hydrothermal waves trav- 
eling in a laterally-heated fluid layer. We have shown 
that, although the regime studied here is rather far from 
the onset of waves, the large-scale modulations of the 
basic pattern are governed by a one-dimensional com- 
plex Ginzburg-Landau equation and we estimated its full 
set of coefficients. This was made possible by showing 
that the localized amplitude holes observed experimen- 
tally correspond to the Bekki-Nozaki hole solutions of 
the CGLE. The overall consistency of our results was 



checked. Since the operating regime of the CGLE at the 
estimated parameter values does not exhibit sustained 
disorder, it would be interesting to analyze experimental 
data collected at other parameter values in the hope of 
reaching spatiotemporal chaos regimes of the type exhib- 
ited by the CGLE. This is left for future work, together 
with an attempt to obtain a better control of the system 
by forcing the behavior of the source. More generally, 
we believe that using the localized structures or defects 
of pattern-forming systems to determine quantitatively 
their relevant amplitude equations can be a rewarding 
approach to this difficult experimental problem. 
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